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ABSTRACT 

The observed spectra of the microlensed sources towards the Galactic bulge may be 
used as a tool for studying the kinematics and extinction effects in the Galactic bulge. 
In this paper, we first investigate the expected distribution of the microlensed sources as 
a function of depth within the Galactic bulge. Our analysis takes a magnitude limited 
microlensing survey into account, and includes the effects of extinction. We show that, 
in the current magnitude limited surveys, the probability that the source lies at the far 
side of the bulge is larger than the probability that the source lies at the near side. We 
then investigate the effects of extinction on the observed spectra of microlensed sources. 
Kurucz model spectra and the observed extinctions towards the Galactic bulge have 
been used to demonstrate that the microlensed sources should clearly show the effects 
of extinction which, in turn, can be used as a statistical measure of the contribution 
of the disk lenses and bulge lenses at different depths. The spectra of the microlensed 
sources provide a unique probe to derive the radial velocities of a sample which lies 
preferentially at the far side of the Galactic bulge. The radial velocities, coupled with 
the microlensing time scales, can thus be useful in studying the 3-dimensional kinematics 
of the Galactic bulge. 

Subject headings: Galaxy: structure — Galaxy: stellar content — stars: kinematics 

1. INTRODUCTION 

The topic of gravitational microlensing has experienced a dramatic increase in scientific interest 
over recent years. This has been largely due to the realization of its wide-ranging applications, such 
as the detection of planets and the study of Galactic structure. More than 400 microlensing events 
have now been detected toward the Galactic bulge and Magellanic Clouds by the microlensing 
survey teams EROS, MACHO, OGLE, DUO, and MOA. The contribution of the bulge stars to the 
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observed microlensing optical depth towards the Galactic bulge was estimated to be approximately 
60%, the other 40% being due to the stars in the Galactic disk (Kiraga & Paczyhski 1994). The 
optical depth to gravitational microlensing, as estimated from the observed microlensing events 
by OGLE (Udalski et al. 1994), MACHO (Alcock et al. 1997), EROS Afonso et al. (2003), and 
MOA Sumi et al. (2003), is significantly larger than the values predicted by theoretical models 
(Paczyhski 1991). It has been argued (Paczyhski et al. 1994) that the observed optical depth can 
be best explained if the effect of the Galactic bar and its inclination are correctly taken into account. 
A consequence of the fact that a large fraction of the events are due to bulge-bulge Icnsing is, as we 
will show in more detail later, that the lensed stars will preferentially be located on the far side of 
the bulge in order that there be sufficient stars along the line of sight to cause microlensing. It was 
suggested by Stanek (1995) that this would mean that there should be a systematic offset in the 
apparent magnitude between observed stars and lensed stars. We present a similar discussion by 
examining the model spectra of different spectral classes and investigating the effects of extinction 
on stars located on the far side of the bulge. 

We first calculate the contributions of the different layers in the bulge to the microlensing 
optical depth, as has previously been calculated by Kiraga &: Paczyhski (1994) and Zhao, Spergel, 
& Rich (1995). We have, however, considered the additional effect of extinction in some detail. 
We then proceed with a discussion on how the spectra of the microlensed sources can be used as a 
measure of the extinction which, in turn, can be used as a statistical measure of the contribution 
of the disk lenses and bulge lenses at different depths. We show that the spectra of the microlensed 
sources can be a useful probe for the 3-dimcnsional structure of the Galactic bulge. We have also 
undertaken such a spectroscopic study of the microlensed sources, the details of which are the 
subject of another paper (Kane k. Sahu 2003). 



2. CALCULATING THE PROBABILITY OF MICROLENSING 

The aim here is to calculate the microlensing probability for sources at various depths within 
the Galactic bulge and to see how the microlensing optical depth varies for sources at various depths 
taking the effect of extinction into account. We first develop the theory for a general distribution of 
microlensing sources and will later restrict the source distribution to those located in the Galactic 
bulge. For these calculations, the formalism used by Sahu (Sahu 1994a,b) shall be adopted. 

Firstly, the number of sources observed at various depths shall be estimated. As we will discuss 
in more detail later, the intrinsic distribution of the Galactic bulge sources can be expressed by 
an exponential distribution, the density being maximum at the central part and falling on either 
side. The extinction and the distance modulus increase as one goes deeper into the bulge, and both 
these effects make the sources fainter. So, in a magnitude limited survey (such as the current ones), 
the combined effect of the extinction and the distance is a monotonous decrease in the observed 
number of stars as we go from the nearest to the farthest region along a given line of sight, which 
needs to be folded in with the intrinsic distribution of the sources in order to obtain their observed 
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distribution. 

To express the above effects, let us assume that the ratio of the stellar number density at the 
central part of the bulge to the stellar density in the near-side of the bulge is ai and the ratio 

of the stellar number density at the central part of the bulge to the stellar density in the far-side 
of the bulge is 02- Thus, ai and 02 include the effects of localized extinctions and the intrinsic 
distributions. 

Let Ntot be the number of stars being monitored in a region. Assuming the extinction to be 
uniform in depth within the bulge, the number of observed stars at any layer dx, at a depth of x 
(as measured from the central plane of the Galactic bulge, see Figure 1), may be expressed as 



-x/d 

Nqo^''''^ ifDs> Rq 



AT r \ ) ^0^1 if < Ro .... 



where No is a constant and represents the observed number of stars per unit depth in the central 
region of the Galactic bulge and d is the depth to which the bulge extends on either side from the 
center. Rq is the distance to the central region of the bulge and Dg is the distance to the source 
(i.e., the distance from the observer to the plane corresponding to x so that x = \Ds — Ro\). 

Integrating this expression along the line of sight, we can write 

rd , fd 



Ntot = I Noa^''''^dx+ I Noa/^'^dx 
Jo Jo 



= No 



In ai \ ai J In 02 V a2 



(2) 



The fraction of area covered by the Einstein rings of all the individual stars lying in front of a 
source at distance Dg may be expressed as 

Af{Ds)= r\RUDa)n{Da)dDa (3) 
Jo 

where is the distance to the lens and n{D^ is the stellar number density at distance D^. If we 
define p{Da) as the average mass density at depth then we see that n^D^) = p{D(i)/M. Since 
the Einstein ring radius for a lens of mass M in this case is given by 



UGM DdiDs-Dd) 



Equation (3) can be rewritten as 



i-Ds 



A,(D.) = -^ p(Di)^^ '-LdDi (5) 
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The instantaneous probability that an observed microlensed star is at a given distance Dg, can 
be written as 

Af{D,)NoUDs) .g. 
Ntot ^ ' 

Note that, in our desire to treat both the source and lens distributions in a similar manner, we 
have chosen our expression for the number of observed stars at Dg to be Nobs{Ds) dDg, and not 
Nohs{Ds) dDg dO,. This would suggest that No{jg{Ds) may increase with Dg since the solid angle 
dO, oc Dg. On the other hand, in a magnitude- limited sample, the Nobs{Dg) would decrease with Dg 
since the observed flux decreases as Dg^. The above two eflFects tend to cancel each other, but not 
exactly, particularly when there is extinction. In our general case of the exponential distribution 
for the bulge (section 2.3), we have chosen ai and 02 to be different, which should absorb this extra 
effect. 

The total microlensing optical depth r can now be calculated. If Ntof stars within the bulge 
are being monitored, then the instantaneous probability of observing a lensing event with an am- 
plification of ^ > 1.34 is given by 

P= Nobg{Dg)Af{Dg)dDg (7) 

Jbn 

where we assume all the monitored sources to be within the Galactic bulge and, for convenience, 
we have defined hn = Ro — d and bf = Ro + d (i.e., the distances from the observer to the near side 
and far side of the bulge respectively). 

The optical depth for microlensing, which is the instantaneous probability that a given star is 
microlensed, can be written as 

T = P/Ntot (8) 

Assuming a uniform extinction with depth is a simple model with some inherent uncertainties 
and is unlikely to be representative of the real extinction in the line of sight to microlensed sources. 
However, this docs not affect the conclusions of this analysis. A full modeling of the extinction has 
its own uncertainties and is beyond the scope of this work. 

It should be noted that in our calculations we are not comparing extinctions at different parts 
in the bulge but rather we are comparing the differential extinctions between microlensed and 
non-microlensed sources in the same part of the sky. It is know that the dust layer when looking 
towards the bulge is patchy and predominantly in the disk (Arp 1965). However, this will not 
significantly affect our calculations as the angular distance between the microlensed source and the 
corresponding lens will be small enough such that the foreground extinction may be considered to 
be the same for both objects (Zhao 2000). 

We will carry out these calculations using 3 different approximations: (i) constant density of 
stars between the observer and the source, (ii) constant (but different) densities for the disk and 
the bulge, and (iii) an exponential density model for the bulge. The last approximation is likely 
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to be the best model, and hence the results of this model will be used in our forthcoming spectral 

analysis (Kane &: Sahu 2003). For completeness however, we present the calculations in the other 
two approximations, since these calculations may be useful for other lines of sight, including the 
lines of sight far from the Galactic Center, and possibly towards the LMC and M31. 



2.1. Constant Density Between Observer and Source 

First, let us assume a constant density of matter for the lenses between the observer and the 
source and make the substitution y = D^/Dg, then equation (5) becomes 



Af{Ds) = ^ p{Dd)y{l-y)dy 



3c2 





(9) 



Since the matter density in the line of sight is constant, extinction is the only effect which 
changes the observed number density of sources in the line of sight. Let us define the intrinsic 
stellar density Nn = Nqu where Nq is the observed stellar density at the mid-plane along the line 
of sight. Let dc = 2d be the total depth, and Uc = which corresponds to the ratio of the stellar 
density from the near to the far side. In this case, equation (1) takes the form 

NobsiPs) = iV„a-^^/'^= (10) 

Integrating this expression yields 



Ntot = 1^ Nr.a-'^^''''^ dDs = Nnd, [^^-^f^ ] (11) 



1 - 1/ac 



Substituting equations (9-11) into equation (6) yields 



3c2d, Vi-IKJ ^ ' 

This probability p is plotted in Figure 2 as a function of Dg assuming a density of p = 0.20 Mqpc"^. 

As seen in Figure 2, for Oc = 1 (which corresponds to no extinction), the probability that 
a given observed microlensed source is at a distance Dg increases monotonically with Dg. It is 
clear that Af{Dg), the fractional area covered by the lenses for a source at distance Dg, increases 
monotonically with distance. However, for Oc > 1, the number of observed stars decreases with 
distance because of extinction. As seen in equation (6), p is a multiplication of these two quantities 
which, as demonstrated in Figure 2, increases steeply with distance for Oc = 1, and increases less 
steeply with distance for increasing values of Qc- This is also demonstrated in Table 1 which gives 
the probabilities that a given microlensed star is at distances Dg = 7.0 kpc, Dg = 8.0 kpc, and 
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Dg = 9.0 kpc for various values of Uc- It is worth reiterating here that p = [Af{Ds)NQhg{Ds)/Ntot] 
is the probabihty that a given observed microlensed source is at a distance Dg, which is distinct 
from the probability that a given observed star at the distance Dg is microlensed; the latter simply 
scales as the microlcnsing optical depth r which monotonically increases with the distance. 

Substituting equation (9) and equation (11) into equation (7) gives 

where the integrand can be solved analytically via repeated integration by parts. Dividing the 
result by N^ot yields an expression for the optical depth for gravitational microlensing 

27rGp4 fle / 2 1 2 2_\ 

Sc^ Qc — 1 \(lnac)^ ttc Oclnttc ac(lnac)^/ 

The value of r as given in equation (14) is plotted in Figure 3 as a function of Oc, which shows 
the variation of optical depth with increasing extinction. As seen from this figure, the optical depth 
in this case is close to 5 x 10"'^ for a large range of ac since the intrinsic matter distribution, which 
plays the major role in the determination of optical depth, is constant. 

However, the assumption of constant density along the line of sight may be valid for microlensed 
sources in the Galactic disk (e.g., OGLE-1999-CAR-Ol (Udalski et al. 2000)), but can not be applied 
to sources in the Galactic bulge. 



2.2. Constant Density for Disk and Bulge 

If we assume a constant density of matter for the lenses in the disk and the bulge, given by pd 
and ph, respectively, equation (5) becomes 

AttG f Dd{Dg-Dd) ,^ , Dd{Dg-Dd) \ 
MDs) = ^[l Pd dD, + j^^ p, dD,^ (15) 

Since the optical depth of microlensing of disk stars is relatively small (Kiraga & Paczyhski 
1994), we have assumed that all the sources arc within the Galactic bulge, i.e., bn < Dg < bf. 

Making the substitution y = D^/Dg, equation (15) simplifies to 

47rGi^2 f ^br./Ds /.I 

'K/Ds 



Af{Dg) = r^ipd yil-y)dy + pj y(l-y)dy (16) 

C \ Jo Jb„/D., / 



Solving this yields the following expression for the fractional area covered by the Einstein rings, 
which is the same as the probability that a given star at a distance Dg is microlensed at any given 
time: 

AttG f (hi hl\ (Dl hi hlW 
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The first term in equation (17) is the contribution of the disk lenses and the second term is that of 
the bulge lenses. 

Since the intrinsic number density of sources within the bulge is constant in this case, extinction 
is the only effect that changes the observed number density of sources in the line of sight. As in 
the previous case, let us define the observed stellar density at the near side Nn = Nqg where Nq is 
the observed stellar density at the mid-plane of the bulge. Let dc = 2d be the total depth of the 
bulge, and Oc = which corresponds to the ratio of the stellar density from the near to the far 
side of the bulge. In this case, equation (1) takes the form 

KM = iVnar (18) 

where Xc (= Ds — bn) is the depth as measured from the near side from the bulge. Integrating this 
expression yields 

Ntot = ^ ' Nna-^^^''^' dXc = Nndc (^^^^) (19) 

Let Ro = 8.0 kpc, 6„ = 7.0 kpc, 4 = 2.0 kpc, pd = 0.07 Mgpc"^, and pi, = 0.50 Mqpc"^. Sub- 
stituting equation (17-19) into equation (6) yields the instantaneous probability that an observed 
star at a given distance Dg is microlensed 

This probability p is plotted in Figure 4 as a function of Ds using the above parameters for the 
bulge. 

Figure 4 shows that, when a constant density is assumed for the bulge, the probability that 
a given microlensed source (in an observed sample of stars) is on the far side of the bulge quickly 
becomes dominated by the level of extinction. Indeed, as shown in Table 2, the probability that 
a given microlensed star is located at the near side of the bulge {Dg = 7.0 kpc) is larger than the 

probability that that a microlensed source is located at the far side of the bulge {Dg = 9.0 kpc) for 
ac ^ 2, simply because too few of the sample of obseerved stars are at 9 kpc. This is analogous to 
the previous case, where the reasons have been explained in detail. 

For completeness, the total microlcnsing optical depth will now be calculated. Substituting 
equations (17-19) into equation (7) yields an expression for the optical depth for gravitational 
microlensing 

-mT^) f: - -'--'"^ (44)- f «-'^-2) ... 

The value of r as given in equation (21) is plotted in Figure 5 as a function of Cc, which shows 
the variation of optical depth with increasing extinction. 
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For Oc = 3.0 the optical depth is ^ 6.76 x 10~ and for Uc = 7.0 the optical depth is 
5.16x10"^. These values are slightly lower than an optical depth of ~ 8.5x10"^ estimated by Kiraga 
& Paczyhski (1994) using a similar Galactic model but without taking extinction into account. This 
is the result one would expect when extinction is added to the optical depth calculations. However, 
our calculation is dependent upon the choice of pd and ph which, given the nature of the model, 
will be crude estimates at best. 



2.3. Exponential Density Model for Bulge 



A more realistic approximation for the density distribution of lenses between the observer and 
the source is to use an exponential- type function to represent the density profile over the Galactic 
bulge which includes the disk component. The model used to do this is a simplified version of the 
E2 triaxial Galactic bulge model fitted by Dwek et al. (1995) with an exponential cutoff at 2.4 kpc. 
This model is defined as 

PE2{i,j,k) = poexp(-r) 



where 



^0/ \3o 



+ — 



1/2 



(22) 



(23) 



is the distance from the Galactic centre. Treating this from a purely radial perspective and centering 
the coordinate system on the observer results in the density profile given by 



p{Dd) = Po exp 



where the normalisation constant is given by 



Po 



\Ro-Dd\ 
io 



M 



S-Kiojoko 



(24) 



(25) 



The value of pQ shall be calculated by adopting the total bulge mass estimated by Dwek et al. 
(1995) of M = 1.3 X 10^° Mq and using the scale lengths of io = 0.71, jo = 0.18, and ko = 0.25 kpc 
for a Galactocentric distance of Rq = 8.0 kpc. This model is valid for a line of sight that passes 
approximately through the Galactic center {l,b) = (0,0). Although it would be interesting to 
conduct a thorough analysis of a line of sight which passes through Baade's window, {I, b) = (1, —4), 
changing the direction slightly will alter the stellar density profile but the extinction model will 
remain largely the same and hence so will the results. 



Using this model, equation (5) becomes 



.4,(0.) = exp 
c Jo 



Jo 



\Ro - Dd\ 
io 



DdjDs - Dd) 



(26) 
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Solving this equation (see Appendix) yields the following expression for the fractional area covered 
by the Einstein rings: 



for Ds < Ro 



where 



(27) 



Following equations (1) and (2), we can express the observed density of sources as 

NobsiDs) = 



Noa2 



{Ds-Ro/d 



if Ds > Ro 



(28) 



where Nq is the the observed number of stars per unit depth at the mid-plane of the bulge. We 
have assumed the lens distribution to be distinct from the source distribution within the bulge since 
the source density distribution has an extra dependence on the sensitivity of observations and the 
extinction. Note that when oi and 02 are greater than unity, the observed source density peaks at 
the mid-plane, which is a realistic case for the Galactic bulge. Also, the definitions of oi and 02 
require that 04 < 02 since the extinction reduces the intrinsic density contrast between the center 
and the near side of the bulge, whereas it increases the density contrast between the center and the 
far side of the bulge. As before, d is the depth to which the bulge extends on either side from the 
center. 

Integrating this expression along the line of sight as in equation (2), we can write the total 
number of observed sources 



Ntot = No 



_d_ A _ ^\ _d_ (^_]_ 
In ai \ ai J In 02 \ 02 



or alternatively 



Ni 



tot 



Nod 



(oi — l)a2 In 02 + (02 — l)ai In ai 
aia2 Inai In 02 



(29) 



Substituting equations (27-29) into equation (6), we can write the instantaneous probability 
that an observed microlensed star is at a given distance Dg 



P 



iirGpo 
4wGpo 



aia2 In ai In 02 



(ai — l)a2 Ina2 + (a2 — 1)11 Inai 

aia2 In ai In a2 

(ai — 1)02 lna2+(o2 — l)ai Inai 



-{Ro-Ds)/d 



-Ro)/d ( io{2ilP-2Rl+Ds-iy 
Ds 



for Ds < Ro 
for Ds > Rq 



(30) 
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where we have assumed that all the sources are within the Galactic bulge, i.e., bn < Dg < bf. The 
probability p is plotted in Figure 6 as a function of Dg using the above model. 

Figure 6 shows that an exponential density model dramatically affects the probability that a 
microlensed source is at a distance Dg. For likely values of ai < 2 and 02 < 2, the probability of a 
microlensed source being in the far-side of the bulge is higher. For higher values of ai and 02, the 
probability peaks at the central plane. This is also seen from Table 3 which gives the probabilities 
that a given microlensed star is at distances Dg = 7.0 kpc, Dg = 8.0 kpc, and Dg = 9.0 kpc for 
ai = 1. The table shows that for this particular case of ai = 1, the probability that a microlensed 
star is on the far side of the bulge is larger than the probability that it is on the near side of the 
bulge if a2 < 3. 

The total microlensing optical depth can now be calculated. Substituting equation (7) and 
equations (28-30) into equation (8) yields an expression for the optical depth for gravitational 
microlensing 



ai02 Inoi Ino- 



(ai - 1)02 In 02 + (a2 - l)ai In oi _ 
Jb„ jRf) Dg 



(31) 



The value of the microlensing optical depth r as given by equation (31) is plotted in Figure 
7 as a function of 02, for various values of oi ranging from 1 to 5. As seen from the figure, 
the microlensing optical depth increases with ai and decreases with a2- For realistic values of 
ai = 2.0 and 02 = 3.0, the optical depth is 4.50 x 10~^, and for ai = 2.0 and 02 = 2.0, the 
optical depth is 4.84 x 10~^. These estimates for the microlensing optical depth towards the 
Galactic bulge are comparable to the values of (3.3 ± 1.2) x 10~^ and (3.9]t]^;2) x 10~^ measured by 
Udalski et al. (1994) and Alcock et al. (1997) respectively. Alcock et al. (2000a) derive a value of 
Ttotal = (2.43lg3g) ^ 10~^ from difference image analysis, and Popowski et al. (2000) derive a value 
of T = (2.43 ± 0.4) x 10^^, which are also consistent with our estimated value. Zhao, Spergel, & 
Rich (1995) and Paczyhski et al. (1994) have estimated a similar value of (2.2 it 0.45) x 10~® and 
(3.3 ± 1.2) X 10^^, respectively, by taking the inclination of the bar into account. 

The limiting apparent magnitude of the current surveys is about my = 20 which at the 
near side of the bulge would correspond to an absolute magnitude of My = 5.8 and an absolute 
magnitude of My = 5.3 at the far side of the bulge were there no extinction. However, the estimated 
foreground extinction towards Baade's window has been estimated by Stanek (1996) to lie in the 
range of 1.26 to 2.79 magnitudes, depending on the line of sight. Hence the limiting absolute 
magnitudes of the stars in this direction lies is in the range My = 4.5 to 2.5 (plus an additional 0.5 
magnitudes because of the larger distance). If the distribution of stars amongst different spectral 
types is assumed to be similar to what is observed in the solar neighbourhood (Wielen, Jahreiss, 
Sz Kriiger 1983; Kroupa, Tout, & Gilmore 1993) then the stellar number density is relatively flat 
in this region. This distribution shows that the ratio in the observed stellar number density from 
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the center to the near or far side of the bulge is between ~ 1.3 to 2.8, depending on the internal 
extinction. 

For oi = 2.0 and 02 = 3.0, the probability that an observed microlensed source is at a distance 
Dg = 9 kpc is ~ 8 times larger than the probability that it is at a distance Dg = 7 kpc. For 
oi = 2.0 and 02 = 2.0, the probability that an observed microlensed source is at a distance Dg = 9 
kpc is also « 13 times larger than the probability that it is at a distance Dg = 7 kpc, and the full 
probability distribution as a function of Dg is shown in Figure 6. This shows that, on average, stars 
that are lensed in Baade's window are more likely to be at the far side of the bulge than the near 
side for reasonable values of ai and a2- 

If lensed stars are indeed predominantly on the far side of the bulge then this provides us 
with a useful tool for studying the Galactic structure at the far side of the bulge. Radial velocity 
measurements from spectra of microlensed sources combined with the measured time scale of the 
events may be used as a unique probe into the 3-dimensional kinematics of the far side of the bulge. 

2.4. Effect of the Galactic Bar 

We have so far neglected the effect of the Galactic bar. The Galactic bar, however, is known 
to exist both from infrared and microlensing observations (Paczyhski et al. 1994; Gerhard 2000; 
Hammersley et al. 2000). A number of photometric and dynamic indications also point out the 
presence of the Galactic bar (Binney et al. 1991; Whitelock, Feast, & Catchpole 1991; Weinberg 
1992). The effect of the bar is not only obvious from the microlensing observations, but its incli- 
nation might be essential to account for the optical depth observed in at least some of the lines 
of sight (also see Binney, Bissantz, & Ortwin (2000)). Furthermore, the inclination of the major 
axis of the bar with respect to the line of sight, as independently determined from the microlensing 
observations, is consistent with the inclination of ~ 15 degrees determined from the earlier infrared 
observations. So a natural question to ask is, how does the presence of the Galactic bar effect the 
analysis presented here? In particular, what is the effect of the Galactic bar on the spectroscopic 
observations, both in terms of extinction and kinematics? 

The effect of the Galactic bar is expected to be seen even more clearly in the spectroscopic 
observations. The reason is two-fold. The first part has to do with the extinction. A good part 
of the Baade's window is thought to have little extinction. The Galactic bar on the other hand, 
which clearly occupies only a part of the region surveyed for microlensing (cf. Blitz & Spergel 
(1991); Udalski et al. (1994); Alcock et al. (1997)), is bright in the far-infrared wavelengths as seen 
in the IRAS SKYFLUX maps (Beichman et al. 1985), clearly indicating the presence of cool dust. 
In such a case, the microlensed sources in the region of the Galactic bar should certainly show 
larger extinction compared to the unlensed sources. The second part has to do with the expected 
kinematics. Although there is some uncertainty in the models, the objects in the Galactic bar, 
in general, are expected to be kinematically distinct from those in the outer parts of the bulge 
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(see, for example, Kent (1987); Blitz & Spergel (1991)). Thus the effect of the bar may be more 

prominently seen in the kinematic distribution. Indeed, the kinematics of the microlensed sources 
may provide very meaningful constraints on the structure of the bar. In order to see this extra 
effect of the Galactic bar, however, the observations must be done for the sources in a restricted 
region covering the Galactic bar. 

2.5. Blending Effects 

The Galactic bulge fields are crowded in general. As a result, many of the microlensed sources 
may be blended with other stars, and this effect must be taken into account in estimating the 
extinctions and the radial velocities of the microlensed sources from the observed spectra. (Note 
that the blending has no effect on the theoretical estimate of the optical depth described above. 
The discussion here refers only to the spectroscopic observations and role of blending in such obser- 
vations). As the recent HST images towards the LMC and the Galactic bulge show, blending can 
be an important effect. In the case of the microlensing events towards the LMC, each 'microlensed 
star' typically splits into 2 or more sources in the high spatial-resolution HST image (Alcock et al. 
2000b). Towards the Galactic bulge, however, for which the distance is about 7 times smaller than 
the LMC, blending may be less severe. 

The effect of blending is discussed in detail by Di Stefano & Esin (1995), including the effects 
of blending when inferring properties of underlying populations through the statistical study of 
Icnsing events (Wozniak & Paczyhski 1997; Dominik & Sahu 2000). The effect of blending can be 
summarized as follows: (i) Blending makes it more difficult to observe a microlensing event since 
the observed amplification is smaller than the actual amplification. This decreases the efficiency of 
the detection of a microlensing event, (ii) Blending enables some stars which are otherwise invisible 
in the sample to be included in the sample of the monitored events. This increases the efficiency of 
detection, (iii) The effect of blending is to increase the number of monitored stars, thus increasing 
the net efficiency of detection. 

The first effect tends to offset the latter two. The net effect depends on the brightness of 
the source/lens, and the crowding of the field. So far as our analysis of extinction is concerned, 
since the blending star is expected to be preferentially closer than the microlensed source, blending 
dilutes the effect of extinction. As Di Stefano k. Esin (1995) point out, the effect of blending is 
more important for fainter sources. A full analysis of blending is beyond the scope of this work. 
However, our preliminary estimate STiggcsts that, if the microlensed star is brighter than niy ~ 18, 
the contribution of the blended star is about 10%. This contribution increases to more than 50% 
for stars with my ~ 20. Thus, if the analysis is confined to brighter sources (as is the case for the 
sources presented in Kane & Sahu (2003)), the effect of blending can be neglected in the analysis 
of the spectroscopic observations. 
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3. 



EFFECTS OF EXTINCTION ON SPECTRA 



It has been shown that a predominant fraction of microlensed sources are located on the far 
side of the Galactic bulge. This can be observationally confirmed by observing the extinction effects 
in the spectra of the lensed stars which, in turn, can be used to estimate the fractional contributions 
of the disk and bulge stars to the total microlensing optical depth. 

To simulate the effects of extinction on various spectral types, the spectral models used were 
those from the Kurucz database (Dr. R. Kurucz, CD-ROM No. 13 (Kurucz 1993)). Each of the 

Kurucz models used have been normalized for a solar metallicity and a distance of Dg = 9.0 kpc. 
Each of the model spectra arc recalculated for Eb^v values of 0.0, 0.2, 0.4, 0.6, 0.8 where the 
extinction corrected spectrum, Sc, is calculated from the raw spectrum, s^, via 



The data for the Galactic extinction, e^, was taken from Seaton (1979). 

Table 4 shows the stellar parameters used for the model spectra, where T^jj is the effective 
temperature, log g is the log gravity, and my is the apparent magnitude of the star. Then varying 
Eb-v from 0.4 to 0.8 is fairly representative of the levels of extinction that exist within Baade's 
window, suggested by Stanek (1996) to be in the range of 1.26 to 2.79 magnitudes. 

Figures 8, 9, 10, and 11 show the effect of extinction on observed spectra of various types of 
stars given in Table 4. There is a significant change in spectral features, such as slope and line 
strength, from an Eb-v value of 0.0 to 0.8. The slope and the features can be used to quantitatively 
estimate the extinction of the microlensed stars in comparison to the general sample in the Galactic 
bulge. It is apparent from the model spectra that, although blue stars in the main sequence will 
always be brighter than red main-sequence stars at all wavelengths, the effects of extinction will 
cause the stellar population at the far side of the bulge to discriminate against blue stars as source 
stars in microlensing. 

From our previous arguments, we expect a majority of the microlensed sources to show ex- 
tinction which would correspond to color excess {Eb-v) values in the range 0.4 to 0.8. Clearly, a 
collection of microlensed source spectra bearing this characteristic would become statistically sig- 
nificant when estimating the contribution of bulge-bulge lensing to the microlensing optical depth. 
If this is shown to be statistically the case then this simple method can be used as a statistical 
distance indicator for microlensed sources. 

A knowledge of the distance to the source enables the distance to the lens to be estimated. 
To demonstrate this, the Galactic bulge model outlined in Section 1.3 shall be used. For a source 
located at a distance Dg, the fractional area covered by the Einstein rings of the intervening stars 



(Sc)A = (5r)AXl0-(^«-^x^«)/2-5 



(32) 



where 



Eb-v 



(33) 
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at a depth dD^ at a distance Dd is 



Af{Dd)dDd 



po exp - 



Ro-Dd\ '] DdjDs - Dd) 



dDd 



(34) 



As discussed earlier, the probability of the source being microlensed is highest when the source 
is at the far side of the bulge. The probability distribution of the location of the lenses (in units 
of the distance to the source) is simply given by Af{Dd), which is shown in Figure 12. The 
figure shows that the probability peaks at a distance of D^, = O.SbDg which, for Dg = 9.0 kpc, 
corresponds to = 7.8 kpc. It is worth noting that this does not imply that most of the lenses 
are at = 7.8 kpc. This only implies that if the source is at Dg = 9 kpc, a lens distance of 0.85 
Dg is twice more likely than a lens distance of 0.6 Dg. 

Since the angular Einstein ring radius is known from the characteristic time scale of the event, 
the extinction exhibited in the spectra of microlensed sources can be used as a means to estimate 
the size of the Einstein ring radius in terms of AU. Detectable deviations in microlensing light 
curves due to planetary masses depend highly upon the star-planet distance in the lensing system 
(Bennett & Rhie 1996; Gaudi &; Sackett 2000), making this method a useful technique to constrain 
calculations of planet detection efficiencies. 



Various Galactic models have been used to demonstrate the theoretical effects of extinction 
upon the microlensing optical depth. The exponential density model provides the best approxima- 
tion of the microlensing optical depth when taking extinction into account. For the exponential 
model, it is shown that stars that are lensed within Baade's window are about 13 times more likely 
to be at the far side of the bulge than the near side. The estimated optical depth is 4.50 x 10~^ 
for oi = 2.0 and a2 = 3.0, which is comparable to the observed optical depths of (3.3 ± 1.2) x lO""^ 
and (3.9+i;2) x 10~^ measured by Udalski et al. (1994) and Alcock et al. (1997) respectively. 

The calculations presented in this paper do not necessarily assume that there is a onoe-to-one 
relation between distance and reddening. Rather we assume that the sources will be more reddened 

than the lens stars and, statistically, other non-microlensed sources in the same narrow field. The 
follow-up study presented in Kane & Sahu (2003), for example, selects comparison stars which are 
very close to the microlensed sources so that the angular distance between them, and hence the 
variation in foreground extinction, is at a minimum. 

The effects of extinction on spectra of microlensed sources have been simulated using Kurucz 
model spectra. It is shown that the majority of microlensed sources should exhibit extinction 
between Eb-v = 0-4 and Eb-v = 0-8. It is also shown that using extinction as a distance 
indicator for microlensed sources may be used to statistically estimate the distance to the lens from 
the probability distribution corresponding to the appropriate Galactic model. 



4. 



CONCLUSIONS 
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The measured extinctions for the microlensed sources can be used to determine the expected 
stellar number density gradient a which, in turn, can be used to determine the optical depth more 
accurately. If there is a smooth velocity gradient within the Galactic bulge, one would expect a 
statistical correlation between the radial velocity and the extinction of microlensed sources which, 
given enough samples, will provide useful information regarding the 3-dimensional velocity structure 
of the far side of the Galactic bulge. Given the variation in extinction through the Galactic bulge, 
the sample would require ~ 100 spectra of microlensed sources in order to clearly show the described 
effects. 

The authors would like to thank Eric Agol for several useful discussions regarding the Galactic 
bulge model calculations. Thanks are also due to Rosanne Di Stcfano for her useful comments 
regarding blending. We would like to thank the referee Hongsheng Zhao for very useful comments 
which led to an improvement of the paper. 



A. FRACTIONAL AREA FOR EXPONENTIAL DENSITY 

This appendix details the evaluation of the integral in equation (26). Let the integral be 
defined as 



/ exp 

Jo 



\Ro - Dd\ 



DdjDs - Dd) 
D, 



dDd 



(Al) 



Due to the absolute value of Rq — D^, this integral is a boundary value problem around the 
Galactocentric distance Rq. Hence, the integral needs to be separated into two integrals defined on 
either of the boundary Rq, as follows 



for Dg < Rq, and 



i: 

L 



exp 



exp 



[Ro-Dd) 



DdjDs - Dd) 
Ds 



dD, 



+ / exp 
JRo 



{Ro - Dd) 
io 

- {D,t - Ro) 



DdjDs - Dd) 
Ds 

D,tjD, - Dg) 
Ds 



dDd 
dDd 



for Ds > Rq. These two integrals will now be considered separately. 



A.l. Source Distance Less Than Rq 



The integral Ii may be expanded into two parts 



Ds 



-Drfcxp 



jRo-Dd) 



Jo Ds 



dDd 



exp 



(i^o - Dd) 



io 



dDd 



(A2) 



Both of these integrals may be evaluated using integration by parts: 



dv f du 

u— ax = uv — / V— ax 
ax / ax 



In the case of the first integral, it is most convenient to define 



u = Dd and 



dv 
dD'd 



exp 



Then it follows that 



■u = exp 

By substituting 7 = (i?o — Dd)/io, v becomes 



= 1 and 



(Ro - Dd 



Dd 



= -io J e '^dj 



io exp 



{Ro - Dd) 



Then using integration by parts 

dv 

u—— dDd = Ddio exp 
dDd 

= ioDd exp 
= io exp 



(Ro - Dd) 

io 

{Ro - Dd) 



lo 

{Ro - Dd) 



lo / exp 



Zq exp 



{Ro-Dd) ' 

io 

{Ro-DdY 

io 



dDn 



{Dd - io) 



In the case of the second integral, it is convenient to define 

{Ro - Dd) 



D, dDd 



exp 



Then it follows that 



du 2Da 



Dd Ds 
Then using integration by parts 



dv 

u——dDd 
dDd 



Dl 
D.. 



iQ exp 



and V = io exp 



(i?o - Dd) 



{Ro - Dd) 
io 



Dd 

io 



^0 

{Ro - Dd) 



2io 
D. 



Dd exp 



{Ro-DdY 



lo 



dDd 



{Ro - Dd) 



2il 



{Ro-Dd) 



{Dd - io) 



{Dl - 2ioDd + 2il) 
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Combining Equation A.3 and Equation A.4 yields 



to exp 



ig exp 



(Rq - Dd) 
io 



Dd-io 



D, 



51 + 



+ Iq exp 



io 



Da=D, 
Dd=0 



D, 



(A5) 



(A6) 



A. 2. Source Distance Greater Than Rq 



Since the first component of integral I2 is identical to integral Ii except for the limits, this 
component may be solved by simply substituting the limits ^ Rq into Equation A. 5. 

Dd=Ro 

inn — #/JI / I J , /.Irs f tJ /'I ^ \ 

h.\ = io exp 



io 



Dj 2ioDd 2i 
Dd — io — + 



■In ■ Ro , 2io-Ro 
zo ito - ^o - + 



2il 

^)-^oexp 



Rq 

io 



io 



2z2 



io{RoD, + 2ioRo + ioe-^/'°{Ds + 2io) - ioD^ - 2il - R^) 



(A7) 



As was the case for integral Ii, the second component of integral I2 may be solved by expanding 
the integral into two parts 



'2.2 



Ds 



Ddexp 



(Dd - Rq) 
io 



dDd 



L 



Ds ^2 
Ro D's 



exp 



(Dd - Rq) 
io 



dDd 



Using integration by parts, the first integral becomes 

(Dd - Ro) 



I2.2.1 

and the second integral becomes 

-^2.2.2 = 



-to exp 



(Dd + io) 



10 



{Dd - Rq 
io 



(Di + 2ioDd + 2i'o) 



so that combining Equation A. 8 and Equation A. 9 yields 
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-^2.2 = io exp 



iQ exp 
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{Ds - Rq) 
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-Dd-io + Tr + 



D 



D. 



+ 



io + 



D, 



Dl 2ioDd , 2il 

Rq ^ 2io-Ro 



D. 
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^o{-Ro-io + ^+ ^ 
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D., 
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Combining Equation A. 7 and Equation A. 10 leads to the expression for integral I2 

/ o / Ds-Rq _Ro \ n ^ Ds-B^o _ Rq 

h = ^ (All) 
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Table 1. Variation of instantaneous probability of microlensing with extinction Oc for a constant 

density between the observer and the source. 





Af{7)NoW^)/Ntot{lQ-^) 


A/(8)iV„6,(8)/iVtot(10-7) 


Af{^)NoU^)/Ntot{lQ-^) 


1.0 


10.92 


14.26 


18.04 


2.0 


8.83 


10.68 


12.51 


3.0 


7.66 


8.85 


9.91 


4.0 


6.87 


7.69 


8.34 


5.0 


6.28 


6.86 


7.26 


6.0 


5.83 


6.24 


6.47 


7.0 


5.46 


5.74 


5.85 
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Table 2. Variation of instantaneous probability of microlensing with extinction Oc for a constant 

density for the disk and bulge. 





Af{7)NoW^)/Ntot{lQ-^) 


A/(8)iV„6,(8)/iVtot(10-7) 


Af{^)NoU^)/Ntot{lQ-^) 


1.0 


3.54 


5.80 


10.22 


2.0 


4.91 


4.02 


3.54 


3.0 


5.83 


3.19 


1.87 


4.0 


6.54 


2.68 


1.18 


5.0 


7.12 


2.33 


0.82 


6.0 


7.61 


2.08 


0.61 


7.0 


8.03 


1.88 


0.47 



- 22 - 



Table 3. Variation of instantaneous probability of microlensing with extinction a2 for an 

exponential density model and ai = 1. 



02 


A/(7)A^<,6,(7)/7Vt„t(10-6) 


Af{%)NoU^)/Ntot{lQ-^) 


A/(9)7V„6,(9)/A^tot(10-6) 


1.0 


0.48 


2.02 


6.07 


2.0 


0.56 


2.35 


3.52 


3.0 


0.60 


2.51 


2.52 


4.0 


0.62 


2.62 


1.97 


5.0 


0.64 


2.70 


1.62 


6.0 


0.65 


2.76 


1.38 


7.0 


0.66 


2.80 


1.20 
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Table 4. Characteristics of example spectral types 



Spectral Type T^ff log 5 my 



MOIII 3800 +1.34 14.4 

GOV 6030 +4.39 19.2 

FOV 7200 +4.34 17.4 

B5III 15000 +3.49 12.4 
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Fig. 1. — Schematic of the geometry used in the Galactic model assuming that all the microlensed 
sources are within the bulge. 
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Fig. 2. — Probability that a given microlensed source is at a distance Dg as a function of Dg, for 
values of ac ranging from Oc = 1 (uniform density of observed stars) to Oc = 10. 
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Fig. 3. — Optical depth as a function of Oc, which is the ratio of the maximum to the minimum 
observed stellar number density in a given hne of sight, assuming a constant matter density for the 
lenses. 
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Fig. 4. — Probability that a given microlensed source is at a distance Dg as a function of Dg, for 
values of ac ranging from Oc = 1 (uniform density of observed stars) to Oc = 10. 
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Fig. 5. — Optical depth as a function of Cc, which is the ratio of the maximum to the minimum 
observed stellar number density in a given line of sight, assuming that the matter density for the 
lenses is constant in the disk and in the bulge. 
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Fig. 6. — Probability that a given microlensed source is at a distance Dg as a function of Dg, for 
values of 02 ranging from 02 = oi to 02 = 10. This has been calculated for (a) ai = 1, (b) oi = 2, 
(c) ai = 5, and (d) ai = 10 with the restriction as expected from physical grounds that ai < a2- 
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Fig. 7. — Optical depth as a function of 02, which is the ratio of the maximum to the minimum 
observed stellar number density in a given line of sight, assuming an exponential matter density for 
the lenses. This has been calculated for ai = 1, ai = 2, and ai = 5 with the restriction as expected 
from physical grounds that ai < 02. 
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Model B5III Spectrum 




4000 6000 8000 10000 

Wavelength (A) 

Fig. 8. — Model spectra of a spectral type B5III star showing the effect of increasing levels of 
extinction corresponding to Eb-v of 0.0 to 0.8 in steps of 0.2. 
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Fig. 9. — Model spectra of a spectral type FOV star showing the effect of increasing levels of 
extinction corresponding to Eb-v of 0.0 to 0.8 in steps of 0.2. 
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Model GOV Spectrum 
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Fig. 10. — Model spectra of a spectral type GOV star showing the effect of increasing levels of 
extinction corresponding to Eb-v of 0.0 to 0.8 in steps of 0.2. 
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Fig. 11. — Model spectra of a spectral type MOIII star showing the effect of increasing levels of 
extinction corresponding to Eb-v of 0.0 to 0.8 in steps of 0.2. 
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Fig. 12. — A probability distribution for lenses given a source located at a distance of Dg = 9.0 kpc. 
In this case, it is most probable that the lens will be located at a distance of = 7.8 kpc. 



